On the Rayleigh-Taylor instability for the incompressible 
viscous magnetohydrodynamic equations 



Fei Jiang'^'*, Song Jiang^, Yanjin Wang*^ 

Institute of Applied Physics and Computational Mathematics, Beijing, 100088, China. 
^School of Mathematical Sciences, Xiamen University, Fujian 361005, China. 



Abstract 

We study the Rayleigh-Taylor instability problem for two incompressible, immiscible, viscous 
magnetohydrodynamic (MHD) flows with zero resistivity and surface tension (or without surface 
tension), evolving with a free interface under presence of a uniform gravitational field. First, 
we reformulate the MHD free boundary problem in a infinite slab as a Navier-Stokes system in 
Lagrangian coordinates with a force term induced by the fluid flow map. Then, we analyze the 
linearized problem around the steady state which describes a denser immiscible fluid lying above 
a light one with a free interface separating the two fluids, and both fluids being in (unstable) 
equilibrium. By studying a family of modified variational problems, we construct smooth (when 
restricted to each fluid domain) solutions to the linearized problem that grow exponentially fast 
in time in Sobolev spaces, thus leading to an global instability result for the linearized problem. 
Finally, using these pathological solutions, we prove the global instability for the corresponding 
nonlinear problem in an appropriate sense. Moreover, we evaluate that the so-called critical 
number indeed is equal to \Jg[Q\/2, and analyze the effect of viscosity and surface tension on the 
instability. 
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1. Introduction 

Considering two completely plane-parallel layers of immiscible fluids, the heavier on top of 
the lighter one and both subject to the earth's gravity. In this case, the equihbrium state is 
unstable to sustain small perturbations or disturbances, and this unstable disturbance will grow 
and lead to a release of potential energy, as the heavier fluid moves down under the (effective) 
gravitational force, and the lighter one is displaced upwards. This phenomenon was first studied 



by Rayleigh [13|, [ij] and then Taylor 16|, and is called therefore the Rayleigh-Taylor (R-T) 
instability. In the last decades, this phenomenon has been extensively investigated from both 
physical and numerical aspects and there are a lot of related results in the literature. In particular, 
many results concerning the linearized problems have been summarized in monographs, see, for 



instance, [2|, Il7[. To our best knowledge, however, there are only few mathematical analysis 
results on nonlinear problems in the literature, due to the fact that in general, passage from 
a linearized instability to a dynamical nonlinear instability for a conservative nonlinear partial 
differential system is rather difficult. 
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The magnetohydrodynamic (MHD) analogue of the R-T instabihty arises when the fluids are 
electrically conducting and a magnetic field is present, and the growth of the instability will be 
influenced by the magnetic field due to the generated electromagnetic induction and the Lorentz 
force. This has been analyzed from the physical point of view in many monographs, see, for 
example, ji], Because of additional difficulties induced by presence of the magnetic field, 
many results concerning the R-T instability of superposed flows could not be directly generalized 
to the case of MHD flows. 

In this paper, we study the R-T problem for two incompressible, immiscible, viscous mag- 
netohydrodynamic (MHD) flows, with zero resistivity and surface tension (or without surface 
tension), evolving with a free interface in the presence of a uniform gravitational field. We will 
prove that in Lagrangian coordinates, the corresponding MHD linearized system is globally un- 
stable in some sense as time increases, and moreover, the original nonlinear problem with or 
without surface tension is globally unstable in some appropriate sense. Next, we formulate our 
problem in details. 

1.1. Formulation in Eulerian coordinates 

We consider the two- fluids free boundary problem for the equations of magnetohydrodynamics 
(MHD) within the infinite slab f2 = x (—1, 1) c M? and for time t > 0. The fluids are separated 
by a moving free interface S(t) that extends to infinity in every horizontal direction. The interface 
divides VL into two time-dependent, disjoint, open subsets VL±{t) so that Vt = f2+(t) UQ-{t) UE(t) 
and S(t) = n+(t) n Cl_{t). The motions of the fluids are driven by the constant gravitational 
field along 63 (the a;3-direction), G = (0, 0, —g) with g > and the Lorentz force induced by the 
magnetic fields. The motion of the fluids is described by their velocity, pressure and magnetic 
field functions, which are given for each t > by, respectively, 

(n±,p±,/i±)(t,-) : fi±(t) -> (M3,M+,M3). 

We assume that at a given time t > 0, these functions have well-defined trace onto S(t). 

The fluids under consideration are incompressible, viscous and of zero resistivity. Hence 
for t > and x = {xi,X2,X3) G ^±{t), the fluids satisfy the following magnetohydrodynamic 
equations: 

dt{g±u±) + div{g±u± ® u±) + dwS± 
divu± = 0, 

dth± + div(u-|- ® h±) — div{h± u±) 
divh± = 0, 

where we have defined the stress tensor consisting of both fluid and magnetic parts by 

\h P 

S± = -fi±{Vu± + Vul) + p±I + ^-^I -h±® h±. 

Hereafter the superscribe T means the transposition and / is the 3x3 identity matrix. The 
positive constants g± denote the densities of the respective fluids. 

For two viscous fluids meeting at a free boundary, the standard assumptions are that the 
velocity is continuous across the interface and that the jump in the normal stress is proportional 
to the mean curvature of the surface multiplied by the normal vector to the surface (cf. ji], [l9|). 
This requires us to enforce the jump conditions 

MIew = 0' (1-2) 

[Su]\j:^,) = kHu, (1.3) 



-9g±e3, 

(1.1) 

0, ^ ^ 
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where we have written the normal vector to as u, /|s(t) for the trace of a quantity / on 
and denoted the interfacial jump by 

[f]\T,{t) ■= /+|s{t) - /-IeW- 

Here we take H to be twice the mean curvature of the surface and the surface tension to be 
a constant k > 0. We will also enforce the condition that the fluid velocity vanishes at the fixed 
boundaries; we implement this via the boundary conditions 

u+{t,x',-l) = 0, u-{t,x',l) = for alH > 0, x' = (xi,X2) G 

Since the fluids are of zero resistivity, the magnetic equations f ll.ip ci are a free transport 
system along the flow, and hence the Dirichlet boundary condition on the velocity at the fixed 
boundary prevents the necessity of prescribing boundary condition on the magnetic field. On 
the other hand, due to the incompressibility f ll.ip /i and also from the physical point of view, we 
assume that the normal component of the magnetic field is continuous across the free interface 
(cf. 0,[l5|) 

[h ■ = 0. (1.4) 

In fact, we will show in the next subsection that the incompressibility of h± and the jump 
condition (11. 4p are satisfied if they hold initially. Therefore, the conditions (ll.ip /i. (II. 4p are 
transformed to the compatibility conditions assumed on the initial magnetic field. 

The motion of the free interface is coupled to the evolution equations for the fluids (II. ip by 
requiring that the surface be advected with the fluids. More precisely, if V{t,x) G denotes 
the normal velocity of the boundary at x G S()f:), then 

V{t, x) = {u{t, x) ■ u{t, x))u{t, x), 

where u{t,x) is the common trace of u±{t,x) onto S(t) and these traces agree because of the 
jump condition (II. 2p . which also implies that there is no possibility of the fluids slipping past 
each other along S(t). 

To complete the statement of the problem, we must specify initial conditions. We give the 
initial interface E(0) = Sq (e.g., Sq = {x G | X3 = 0}), which yields the open sets ^^±(0) on 
which we specify the initial data for the velocity and magnetic field 

(M±,/^±)(0,■):^^±(0)^(M^M3). 

Thus the initial datum of the pressure po can be defined by g±, Sq, u±{0,-) and /i±(0, ■). 
To simply the equations we introduce the indicator function x and denote 

= Q+Xn+ + Q-Xn-, u = u+xn+ + U-Xn^, 

h = h+xn+ + h_xn-, P = P+Xn+ + P-Xn., 

and also define the modified pressure by 

\h\^ 

P = P+ — + 90x3. 

Hence, the equations (II. ip can be rewritten as 

{gdfU + Qu ■ Vu + Vp = h ■ Vh + /iAm, 
dth + u- Vh-h-Vu = Q, 
divM = div/i = 

in 0\E(t) for each t > 0, and the jump condition (II. 3p becomes, setting [g] = g+ — g-, 
[{pi - ix{Vu + Vu^))i']\ = g[g]x3U + h ■ u[h]\„, + kHv. 
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1.2. Formulation in Lagrangian coordinates 

The movement of the free interface and the subsequent change of the domains ^l±{t) in 
Eulerian coordinates will result in severe mathematical difficulties. To circumvent such difficul- 
ties, we switch our analysis to Lagrangian coordinates, so that the interface and the domains 
stay fixed in time. To this end, we define the fixed Lagrangian domains = x (0, 1) and 
n_ = X (— 1, 0), and assume that there exist invertible mappings 

7]l:Q±^ n±(0), 

such that So = ri^{{x3 = 0}), {x^ = 1} = ri'^{{xs = 1}) and {x^ = —1} = ({.X3 = — 1}). 
The first condition means that Eq is parameterized by the mapping rj^ restricted to {x^ = 0}, 
while the latter two conditions mean that 77° map the fixed upper and lower boundaries into 
themselves. Define the flow maps r)± as the solution to 

j dtri±{t,x) =u±{t,ri±{t,x)) , . 

We denote the Eulerian coordinates by (t, y) with y = ri{t, x), whereas the fixed {t, x) G M"*" x Q 
stand for the Lagrangian coordinates. In order to switch back and forth from Lagrangian to 
Eulerian coordinates, we assume that ri±{t, •) are invertible and fl±{t) = ri±(t, fl±), and since u± 
and Tjj. are all continuous across {^3 = 0}, we have E(t) = r]±{t, {x^ = 0}). In other words, the 
Eulerian domains of upper and lower fluids are the image of ft± under the mappings r)±, and the 
free interface is the image of {xs = 0} under the mappings r]±{t, ■). 

Setting r] — X+V+ + X-V-i we define the Lagrangian unknowns by 

{v,q,b){t,x) = {u,p,h){t,r]{t,x)), {t,x) e x Q. 

Defining the matrix A := {Aij)sx3 via = [Dr])^^ := {d.jr]i)^^^, and the identity matrix 
/ = (/ij)3x3, thus in Lagrangian coordinates the evolution equations for rj, v, q, b read as (writing 
dj = d/d^.) 

' dtfji = Vi 

gdtVi + AjkdkTij = bjAjkdkbi, 
< AjkdkVj = 0, (1.6) 

dtbi = bjAjkdkVi, 
< Ajkdkbj — 0, 

where the stress tensor of fluid part in Lagrangian coordinates, T{v, q), is given by 

Tij = qlij - ii{AjkdkVi + AikdkVj). 

Here we have written I^j for the ij-component of the identity matrix I and used the Einstein 
convention of summing over repeated indices. 

To write the jump conditions, for a quantity / = /±, we define the interfacial jump by 

[/l '■— f+\x3=0 ~ f-\x3=0- 

Then the jump conditions in Lagrangian coordinates are 

M = 0, Ibjrijj = 0, {TijUj} = g[g]mni + MbjUj + kHth, (1.7) 
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where we have written n := (ni, ?t,2, "^s) = i-e. 

diT] X d2ri Aes 



n 



\dir] X d2r]\ \Ae3 



{X3=0} 



for the normal vector to the surface = r]{t, {x^ = 0}) and H for twice the mean curvature 
of Since is parameterized by rj, we may employ the standard formula for the mean 

curvature of a parameterized surface to write 

^ _ / \diri\^d^r] - 2{dir] ■ 8211)81827] + \d2r]\^dfr] \ _ ^ 
V |9i?7P|9277|2 - |(9i?7 ■ c?2?7|2 J 

Finally, we require the no-slip boundary conditions 

v-{t,x',-l) = 0, v+{t,x',l) = for alH > 0, x' G (1.10) 

1 . 3. Reformulation 

In this subsection we reformulate the free boundary problem (11. 6p . (I1.7P and (ll.lOp . Our goal 
is to eliminate b by expressing it in terms of rj, and this can be achieved in the same manner as 
in [isj. For the reader's convenience, we give the derivation here. 

Applying An to (ll.6P /i. we have 

Aii8tbi = bjAjk8kViAii = bjAjkdt{dkrji)Au = -bjAjk8kr]idtAii = -bi8tAii, 
which implies that 8t{Ajibj) = 0. Hence, 

AA = 4&;, (1.11) 

k = 8ir],A%bl (1.12) 

Hereafter, the superscript means the initial value. 

With help of (I1.12p . and the geometric identities J = J°, 8k{JAik) = where J = \Dri\, we 
can apply Aik8k to fll.l2p to deduce that the divergence of b (i.e. (11. 6p ^) satisfies 

Ak8kbi = j^A,A{8ir],A'jfi'j) = j^8k{JA,k8iv,A'jib^^) = ^8u{J''A%b]) = A],8kbl (1.13) 

Next, we evaluate the jump It is easy, recalling (II. Sp . to verify that Ae^ is continuous 

across the free interface. Hence, one has 

[6,n,l = l8iV,AlblA,,l^^ = lAlAij^^ = I^^S^- (l-^^) 

Thus, if one assumes the compatibility conditions on the initial data 

A%8kb] = ^, [&>;i=0, (1.15) 
then from (I1.13P and (I1.14p one gets 

Ajk8kbj = 0, {bjUjl = 0. (1.16) 
Moreover, for simplifying notation, we assume that 

^m/^m — with B being a constant vector. (1-17) 
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We remark that the class of the pairs of the data 17°, iP that satisfy the constraints fll.lSp and 
f ll.lTp is quite large. For example, we chose r]*^ =Id and IP =constant vector, then 77° and IP satisfy 
( ll.lSp and fll.l7p . consequently, the pair (77, h) which is transported by the flow will satisfy ( I1.16p . 
Now, in view of (11. lip . (I1.12p . (I1.17p . we represent the Lorentz force term by 

h.A^M, = dir^,AlihlA,A{drV^Alhl) = Al,bl^,i^rV^A',X) = BiB^^^r^,. 

Hence the equations (II. 6p become a Navier-Stokes system with the force term induced by the 
flow map rj: 

{dtVi = Vi 
gdtv^ + AjkdkT.j = BiBrndl^T],, (1.18) 
AjkdkVj = 0, 

where the magnetic number B can be regarded as a vector parameter. Accordingly, the jump 
conditions (II. 7p become 

1^1 = 0, {TijUj} = glglvsni + BiBmldiriijdm'njnj + nHui. (1.19) 

Note that we implicitly admit that BmdmVj^j is continuous across {x^ = 0}. In fact, this follows 
from the assumptions (ll.lSp . (I1.17p . 

Finally, we require the boundary conditions (ll.lUp . 

1.4. Linearization around the steady state 

The system (I1.18p . (ll.lQp and (ll.lOp admits the steady solution with v = 0, rj = Id, 

q =constant with the interface given by //({xs = 0}) = {0:3 = 0}, and hence, n = 63 and 

A = I. Now we Inearize the equations (I1.18P around such a steady-state solution, the resulting 
linearized equations are 

{dtV = v, 

gdtv + Vq = fiAv + BiB^df^r], (1.20) 
divf = 0. 

The corresponding linearized jump conditions are 

M = 0, Iql - /i(Vt; + Vt;^)]e3 = {9[g]v3 + ^^A,'V3)e3 + B^Bildivl (1.21) 

while the boundary conditions are 

v.{t,x',-l) = 0, v+{t,x',l) =0. (1.22) 

Now, we briefly review some of the previous results on the nonlinear Rayleigh- Taylor insta- 
bility without the magnetic fleld. In 1987, Erban jij proved the ill-posedness of the equations of 
motion for a perfect fluid with free boundary. Then, he adapted the approach of jsf to obtain the 
nonlinear ill-posedness of both R-T and Helmholtz instability problems for two-dimensional in- 
compressible, immiscible, inviscid fluids without surface tension 0]. In 2011, for two-compressible 
immiscible fluids evolving with a free interface, Guo and Tice made use of flow maps (cf. (II. 5p ) 
to transfer the free boundary into a flxed boundary and established a variational framework 
for the nonlinear instability in joj, in which with the help of the method of Fourier synthesis, 
they constructed solutions that grow arbitrarily fast in time in Sobolev spaces, leading to the 
ill-posedness of the perturbed problem in Lagrangian coordinates. This is in some sense a com- 
pressible analogue to the ill-posedness of the R-T instability problem for incompressible fluids 
given in jsj. Unfortunately, the approaches in both 0] and ^ could not be directly applied to 
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the viscous flow case, since the viscosity can bring some technical difficulties in the study of the 
nonlinear R-T instability. Hence, Guo and Tice only investigated the stabilized effect of viscosity 
and surface tension to the linear R-T instability (see 0]) for compressible flows, and the corre- 



sponding nonlinear instability still remains open. Recently, Jiang, Jiang and Wang [10[ showed 
the nonlinear instability in some sense for two incompressible immiscible fluids with or without 
surface tension in Eulerian coordinates. 

Unfortunately, the results concerning the nonlinear R-T instability in 0, [g], [io[ are still not 
generalized to the case of MHD flows, because of additional difficulties induced by presence of the 
magnetic field. For two incompressible immiscible fluids evolving with a free interface (the density 
is discontinuous across the interface) under presence of a magnetic field, it was ffist shown by 
Kruskal and Schwarzschild [IH that a horizontal magnetic field has no effect on the development 
of the R-T instability for the linearized equations in the whole space in Eulerian coordinates. 



Recently, for the case of finite slab domain without surface tension, Wang [18| obtained the 
critical magnetic number (denoted by \B\c) for the linear R-T instability. Namely, he gave an 
instability criterion to the linearized problem fll.20p in terms of the value of the magnetic number, 
and pointed out the linear R-T instability in the case B = (0, 0, B) with \B\ < \B\c. In particular, 
unlike 2D case, he also showed that the linearized problem is unstable for any initially horizonal 



magnetic field B = {B, 0, 0). Form [18[, we see that the stabilized effect of the magnetic field is 
more remarkable than that of the surface tension k which only stabilizes the frequency interval 
(0) a/sM7^) for S'liy «^ > 0, also see 0]. Adopting and modifying the approaches in both 
and [18|] , Duan, Jiang and Jiang recently showed the ill-posedness of the linearized system fll.20p 
with fj, = and B = (5,0,0) in the sense of Hadamard and moreover, the ill-posedness of 
the corresponding nonlinear problem in some sense. Finally, we mention that Hwang derived 
the nonlinear instability around different steady states for both incompressible and compressible 
inviscid MHD flows with continuous density, thus extending the results in Q to the case of MHD 
flows with continuous density. 

In this paper, we will study the R-T instability for the problem f ll.lSp . ( I1.19P and f ll.lOp . 
We will prove that the corresponding linearized system f ll.20p -( n..22p with any B = (5,0,0) 
is globally unstable. Moreover, the original nonlinear problem f ll.lSp with or without surface 
tension is globally unstable in some sense. For this purpose, we assume that k > and that the 
upper fluid is heavier than the lower fluid, i.e., 

g-i-> ^ [g] > 0. 



In addition, we will compute that the so-called critical number in 18 1 indeed is equal ^^/g[Q\/2 
(see Remark 13.11 for the details), and analyze the effect of viscosity and surface tension on the 
linear instability (see Remark 13. 2p . The results of the current paper extend the ones in (sl. isl . 



The rest of this paper is organized as follows. In Section 2 we state our results concerning 
the linearized system fll.20p and nonlinear system fll.lSp . i.e. Theorems 12.11 and 12.21 respectively. 
In Section 3 we construct the growing solutions to the linearized equations, while in Section 4 
we prove the uniqueness of the linearized problem and Theorem 12. 1[ In Section 5, we prove the 
global instability of order k to the corresponding nonlinear problem, i.e.. Theorem 12. 2[ 

2. Main results 

Before stating the main results, we introduce the notation that will be used throughout the 
paper. For a function / e L^(fi), we define the horizontal Fourier transform via 



fX^,xs)= [ f{x',x,)e-'''<dx', 

JR2 
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where x',^ G and "■" denotes scalar product. By the Fubini and Parseval theorems, we have 



\fix)\'dx 



47r2 



d^dxs. 



We now define a function space suitable for our analysis of two disjoint fiuids. For a function 
/ defined on Q, we write /+ for the restriction to = x (0, 1) and /_ for the restriction to 
Q_ = X (—1, 0). For s G M, define the piecewise Sobolev space of order s by 



(2.1] 



endowed with the norm 
be 



2 _ 
k 

i=o 

k 

i=o 



+ 



)• For A; G N we can take the H*-norm to 



2\k-j 



'xl± 



2\k-j 



2 



d,^da;3 



for /_ = (—1, 0) and /+ = (0, 1). The main difference between the piecewise Sobolev space H*(f2) 
and the usual Sobolev space H'^{Q) is that we do not require functions in the piecewise Sobolev 
space to have weak derivatives across the interface {x^ = 0}. 

Now, we are in a position to state our first result, i.e. the result of global instability for the 
linearized problem f ll.20p . 



Theorem 2.1. Assume the constant vector B = (5,0,0), then the linearized problem M.2U\) with 
the corresponding jump and boundary conditions M.21\) . is globally unstable in &{Q) for 

every k. More precisely, there exists a constant Ci > 0, and for any k, j & N with j > k and 
for any a > 0, there exist a constant Cj^k depending on j and k, and a sequence of solutions 
{{Vn,Vn,(ln)}'^=i to U.20\) Satisfying the corresponding jump and boundary conditions U.21]) . 
n^) . such that 

1 



hn(0)||H. + ||t^n(0)||H. + |kn(0)||H. < 



n 



but 



ll'^ri(^) llflfe ^ cfc for all t >tn '■= Cj^k + Ciln(an). 
Moreover, there exists a positive constant Aq such that 

ll'7n.('^) IIh'= ^ '^o||'^n(^)||H'= '^^^ ll'^n(^) ||h'= " ^ OO flS t — )■ OO. 



(2.2) 
(2.3) 

(2.4) 



Theorem 12.11 shows globally discontinuous dependence of solutions upon initial data. The 
proof of Theorem 12.11 is motivated by 0] under nontrivial modifications, and its basic idea is 
the following. First, noticing that the coefficients of the linearized equations depend only on the 
vertical variable X3 G (—1,1), we seek "normal mode" solutions by taking the horizontal Fourier 
transform of the equations and assuming that the solutions grow exponentially in time with the 
factor e'*'^'^^''''^^'-'*, where ^ G is the horizontal spatial frequency and A(|^i|, |^2|) > 0. This 
reduces the equations to a system of ordinary differential equations (ODE) with A(|^i|, |^2|) > 
for each ^. Thus, solving the ODE system by the modified variational method, we show that 
A( 1^1 1, 1^2!) > is a continuous function on some symmetric open domain A C G | |^| G 
(0, |^|c)} (see Proposition 13.71) . the normal modes with spatial frequency grow in time, providing 
a mechanism for the global R-T instability, where = ^/qIqI/k, if k > 0, otherwise |^|c = 00. 
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In particular, we can restrict ^ in some symmetric closed sector domain D, such that A(|^|) has a 
uniform lower bound for any ^ G D (see Remark 13.31) . Then, we form a Fourier synthesis of the 
normal mode solutions constructed for each spatial frequency to give solutions of the linearized 
equations that grow in time, when measured in for any k > 0. Finally, we exploit the 

property of the trace theorem to show a uniqueness result of the linearized problem (i.e. Theorem 
14.11) . Our results show that in spite of the uniqueness, the linearized problem is still globally 
unstable in for any k. 

With the linear global instability established, we can show the global instability of the cor- 
responding nonlinear problem in some sense (i.e. Theorem 12. 2p . Recalling that the steady state 
solution to fll.20p is given by ?7 = 0, f = 0, g =constant, we now rewrite the nonlinear equations 
(ll.lSp in the form of perturbation around the steady state. Let 

1] = ld + f] is invertible, q = constant + q, v = + v, A = I — G, 

where 

= I - {I + Dfi)-\ 
Then, the evolution equations fll.lSp can be rewritten for fj, v, q as 

dtfj = v, 

gdtVi + {Ijk - Gjk)dkfij = BiB^df^fji, i = l,2, 3, (2.5) 
divv - ti{GDv) = 0, 

where tr(-) denotes the matrix trace and 

fij = qlij - fi{{Ijk - Gjk)dkVi + {lik - Gik)dkVj). 

The jump conditions across the interface are 

{v} = 0, {fijUj} = g[g]fi-ini + BiB„,ldifiij{Iij + dmVj)nj + nHui, (2.6) 

where n and H are respectively given by (ll.Sp and (11. 9p with rj = ld + fj. Finally, we require the 
boundary conditions 

v-{t,x',-l) = 0, v+{t,x',l) = 0. (2.7) 

We collectively refer to the evolution, jump, and boundary equations (12. 50 - 02. 7p as "the 
perturbed problem" . To shorten notation, for k > we define 

\\iV,v,q){t)U>. = ||r/(t)||H. + ||^^(t)||H^ + Wmih^- 

Definition 2.1. We say that the perturbed problem has global stability of order k for some k > 3 
if there exist 6, C2 > and a function F : [0, 6) satisfying F{z) < C2Z for z G [0, 5), so 

that the followings hold. For any T > 0, fjo, and Vq satisfying 

\\{Vo,vo)\\fik < 6, 

there exist fi,v,q E L°°{0, T; H^(f2)), so that 

(1) {fi,vm = {fio,vo); 

(2) rj, V, q solve the perturbed problem on (0,T) x Q; 

(3) it holds that 

sup \\{fi,v,q){t)\\ii3 < F{\\{fio,Vo)\\nk); 

0<t<T 

(4) r]j- G C'^{Q±), respectively, when k > 0. 

Remark 2.1. In the above definition, in view of (II. 8p and (11.90 . we have admitted that 
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n = 



dir]± X d2i]± 

\diV± X d2V±\ {,3=0} 



and 



H 



dir]±\'^d^r]± - 2{dir]± ■ d2r]±)did2ri± + |^2^±P^f^± 
\dir]±\^\d2r]±\'^ - \dir]± ■ d2T]±\^ 



) 



■ n 



{^•3=0} 



with ri± = Id + r^-i-. For convenience in the subsequent proof of Theorem \2.2\ we will still use rj 
to denote the both //_ and rj+ at the interface {0:3 = 0}. Similarly, fj in fl2.6p also includes the 
both cases of fj^ and fj^ at {xs = 0} except for the notation |(9/f/j]. We point out that Theorem 
12.21 below still holds if fj only denotes fj+ (or fj^) at {x^ = 0}. 

We can show that the property of global stability of order k cannot hold for any k > 3, i.e., 
the following Theorem 12.21 which will be proved in Section 5. 

Theorem 2.2. Assume the constant vector B = {B,0,0), the perturbed problem does not have 
property of global stability of order k for any k > 3. 

The basic idea in the proof of Theorem 12.21 is to show, by utilizing the Lipschitz structure 
of F, that the global stability of order k would give rise to certain estimates of solutions to the 
linearized equations (11.181) that cannot hold in general because of Theorem 12. 1[ We will adapt 
and modify the arguments in to prove Theorem l2.2[ Compared with the perturbed problem in 
[ol. Theorem 5.2], the main difficulty lies in the convergence of the jump conditions of the rescaled 
pressure function sequence, because we could not obtain the strong convergence of the rescaled 
pressure as in To circumvent such difficulty, similarly to fio'], we apply the techniques of 
integrating by parts to avoid using the strong convergence of the rescaled pressure. Also, this 
idea is used in the proof of the uniqueness of solutions to the linearized equations fll.20p in Section 
4. 

Remark 2.2. Theorems 2.1 and 2.2 show that a horizontal magnetic field can not prevent the 
linear and nonlinear R-T global instability in the sense described in Theorems 2.1 and 2.2. From 
the proof we easily see that Theorems 12.11 and 12.21 still hold with B = (0, B, 0) in place of 
B = (5,0,0). In addition, based on the results in we can show that Theorems 12.11 and 
12.21 still remain valid with B = (0,0,5) in place of B = (5,0,0), if one of the following two 
conditions holds: 

(1) K = and \B\ < g[g]/2; 

(2) K, > and |5| sufficiently small. 

3. Construction of a growing solution to linearized equations 

3.1. Growing mode ansatz and the horizontal Fourier transform 

We wish to construct a solution to the linearized equations fll.20p that has a growing H'' 
norm for any k. We will construct such solutions via Fourier synthesis by first constructing a 
growing mode for fixed spatial frequency. 

To begin, we assume a growing mode ansatz (cf. [2]), i.e.. 



v{t,x) = w{x)e \ qit,x) 



■w{x)e *, Tjit.x) = <;(x)e' 



xt 
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for some A > 0, where w = {wi,W2,W3). Substituting this ansatz into fll.20p . ehminating q by 
using the first equation, we arrive at the time-independent system for w and w. 

\gw + Vw = ^Aw + X'^BiBmdf^w, , 
divw = 0, \ ■ ■> 

with the corresponding jump conditions 

M = 0, ItuI - fi{Dw + Dw^)je, = X-\g[g]ws + KA^'Ws)e, + X-'BsBildiwl 

and the boundary conditions 

w{t,x',-l) = 0, w{t,x',l) = 0. 

We take the horizontal Fourier transform of Wi, W2, in fl3.ip . which we denote by either " 
or J^, and fix a spatial frequency ^ = (^i, ,^2) ^ I^^- Define the new unknowns (p^x^) = iwi{^, X3), 
9{x3) = iw2{^,X3), ip{x3) = t&3(^,X3) and 7r(x3) = W3{^,X3), so that 

J^(divw) = ^iV^ + 6^ + 

where ' = d/dx^. Since we consider the case B = {Bi, B2, B3) = {B,0,0), then for if, 9, ip and 
A = A(^) we arrive at the following system of ODEs. 



^ X'g^ - A^ivr + A/x(|e| V - v") + B^Hv = 0, 
X^g9 - A^TT + Xfi{\^\^9 - 9") + B^eiO = 0, 
X^Qi^ + Avr' + A/^(|eP^ - O + B^eii^ = 0, 



(3.2) 



along with the jump conditions 

M = m = m = 0, 

[A/x(ei^ - v^Ol = IA/i(6^ - = 0, (3.3) 
l-2Xfif + X7fj = {g[g]-km^, 

and the boundary conditions 

^(-1) = ^(1) = 9{-l) = 9{1) = V'(-l) = V(l) = 0. (3.4) 

Eliminating tt from the third equation in (13. 2p . we obtain the following ODE for 

- x'pm'^p - r) = A/x(iei V - 2iei v + r') + B^.m^^p - r) (3.5) 

along with the jump conditions 

m = m = 0, (3.6) 

[AMieP^ + V^")l=0, (3.7) 

ixfiir - 3iei V)i = ix'Qn + {9[p] - m'm'i^, m 

and the boundary conditions 

V^(-l) = V^(l) = V^'(-l) = V^'(l) = 0. (3.9) 
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Before constructing a growing solution to ODE of f l3.5p - fl3.9p in next subsection, we should 
introduce so-called admissible set for a growing solution. 
First we define the maximal frequency by 



I \/ q\q\I if K > 0, , , 

\i\c:={ l'^'^' ; (3.10) 

-|-oo, it K = U, 

and the generalized critical magnetic number (depending on the frequency) by 



Bo(lfl):= =up ^M-Aimm fo,|^|<|j|„ (3.U) 



and the generalized critical frequency function by 



•5(161, 161) 



sup 1 3.12 

/^(l6l,l6l)/il^|Mi/ 



for < lei, and < /(O < 5c(|6), where f{\U 161) = I6i?l/I6 and S{\U 161) := oo when 
= 

By [18|, Lemma 3.2], we see that the supremums in both (13.111) and (I3.12p are achieved for 
each ^ and B with ^iB ^ 0. So, (13. lip and (I3.12p make sense. Moreover ^d^il, |,^2|) — > oo as 

/(I61, 161)^0. 

Now, we define an admissible set for a growing solution 

A = G I 16 < I6c, < /(l^il, 161) < 5c(l6), < 16 < ^(161, 161)}- (3.i3) 

By construction, the admissible set possesses the following properties, which will be useful in 
Subsections 3.2-3.4. 

Proposition 3.1. Assume A is defined by Ii3.13\) . then 

(1) the set A is a open set in M^; 

(2) the set A is symmetric on x-axis and y-axis in , respectively; 

(3) there exist three positive constants di < d2 and d^, such that 

D := G I rfi < 16 < d2, l^il < 4} C A. (3.14) 

Proof. The proof is very elementary and we omit it here. 

In the next subsection, we will show that for any fixed B and any ^ G A, there is a nontrivial 
solution ip with A > to the problem (l3.5p -f lX^ . The next proposition shows that -Bc(l6) can 
be in fact represented hy g, [g], k and |6- 



Proposition 3.2. Let 5^(16) be defined by [3ll\) . Then, 5^(16) = ^/MoF^^AW)!^ ■ 
Proof. By virtue of (13. lip , it suffices to prove that 
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For any ip G Hq{—1, 1), we have ip{0) = J^^ip'dy = tjj'dy, which gives 



2 

Hence, by the Holder inequahty, we see that 



1^(0)1 <^ f'mdy. 



< 1 < i. (3.16) 



On the other hand, it is easy to check that the function 



1 + x, X e (-1,0], 
1 - X, X e (0, 1) 

satisfies ips ^ Hq{—1, 1). Furthermore, a simple computation results in 

which combined with f l3.16p implies f lS.lSp . □ 

Remark 3.1. We mention that the definition of critical number and critical frequency (both are 
independent of frequency ^) was introduced by Wang 18j, where he defined the critical number 

by 

\B\^:= sup 



which is equal to ^/g[g]/2 by virtue of Proposition 13.21 

3.2. Construction of a growing solution to ODE 

Throughout this subsection we assume that ^ G A. Similarly to 0, [lH, we can apply the 
variational method to construct a growing solution of fl3.5p - fl3.9p for given ^ E A. For the reader's 
convenience, we sketch the procedure of the construction. 

First of all, in order to apply the variational method, we remove the linear dependence on A 
in fl3.5p - fl3.8p by defining the modified viscosities fl = s/i, where s > is an arbitrary parameter. 
Thus we obtain a family (s > 0) of modified problems 

- x'pim - r) = ^/i(iei V - 2\er + ri + B^^iim - r) (3.17) 

along with the jump conditions 

M = m = 0, (3.18) 

Mlel'^ + ^")1 = 0, (3.19) 
ls^^{r' - 3|e| V)l = IX'gn + {9[p] - ^l^mi'^P, (3.20) 

and the boundary conditions 

V^(-i) = v^(i) = v^'(-i) = V^'(l) = 0. (3.21) 
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Notice that for any fixed s > and C,, f l3.17p -( !3.2ip is a standard eigenvalue problem for — A^. 
This allows us to use the variational method to construct solutions. We define the energies by 

^-^-1 (3.22) 



and 



^^ = 2/ /(l^lV+lV^'Dcixs (3.23) 
which are well-defined on the space Hl{—1^ 1). We define the admissible set for the energy fl3.22p 

A = {i^eHl{-l,l) I J(^) = l}. 
Thus we can find the smallest — by minimizing 

-A2(|^i|,|6|) = «(|6|,|6|):= inf W). (3.24) 

In fact, we can show that a minimizer of (13.240 exists and satisfies the Euler-Lagrange equations 
equivalent to flXT7D - flX^ . 

Proposition 3.3. For any fixed ^ E A, E achieves its infinimum on A. In addition, let ip be 
a minimizer and — A^ := E{%1)), then the pair (ip, \^) satisfies (3. 17) along with the jump and 
boundary conditions /I3.18\) - f3.21]) . Moreover, ip is smooth when restricted to (—1,0) or (0,1). 



Proof. We can follow the same proof procedure as in 18|, Proposition 3.1] to obtain Proposition 
13. 1[ Hence, we omit the details of the proof here. □ 

In order to emphasize the dependence on s G (0, 00), we will sometimes write 

H\um = x{\u\us) 

and 

a{\Um=a{s). 

Next, we want to prove that there is a fixed point such that X = s. To this end, we first give 
some properties of a{s) as a function of s > for each fixed ^ G A. 

Proposition 3.4. Let a{s) : (0, 00) — > M 6e given by ^3.24\ ), then the following conclusions hold: 

(1) a{s) G C[°'^(0, 00) is strictly increasing; 

(2) there exist two positive constants ci := ci{g, g±) and C2 := C2{\i\,f^±), such that 

a{s) > -cil^l + SC2; 

(3) for each fixed ^ G A, there exist two positive constants C3 := C3{\B\,\^\, g, g±, n) and 
C4 := C4(|i?|, \^\,g, g±, /i±, k), such that 

a{s) < — C3 + SC4. (3.25) 
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Proof. The first two assertions can be shown in the same way as in 18, Lemma 3.5] and 0, 
Propostition 3.6]. It remains to prove the third assertion to complete the proof. 
First, observe that the energy E{iIj) can be decomposed in the form: 



where 



E{ij) = \^\'Eo{ij) + sE,{ij), 



1 1 



Then, since ^ G A by virtue of the definitions fl3.12p and fl3.13p . there exists ip G Hq such that 

+ m')dxs < {g[0] - V^^(O). (3.26) 

On the other hand, C^{—1,1) is dense in Hq{—1,1), then there exists a function sequence 
V'n G (-1, 1), so that 

i^n^i' stronly in H^{-1, 1), (3.27) 
which, together with the compact embedding theorem, yields 

?/'„^V^ stronly in C0(-1,1). (3.28) 

Putting fl3.26p -( l3.28p together, we see that there exist a subsequence ipno G Hq{—1, 1), such that 
ipno ^ and 

f (161, mm'l^nol' + Kf)dx, < {g[g] - k\^\') <(0), (3.29) 



which implies 
Thus, we have 



^o(^no) < 0. 

a{s) = inf E{tlj) = inf ^ 

^ Wno) ^ |^|2 -^o(V^«o) ^ „ -^l(^no) ^ 
~J(V^„J ^(^no) ^(V^no) ' ^ ^ 



(3.30) 



for two positive constants C3 := C3(|i?|, |^|, (7, f?±, k) and C4 := C4(|i?|, \'i\,g, g±, fi±, n). □ 

Given ^ G A, by virtue of f l3.25p . there exists a sq > depending on the quantities \B\, \C,\,g, g±, 
fi± and K, so that for s < sq, a{s) < 0. Therefore, we can define the open set 

S = a'\-oo,0) C (0,00). 

Note that S is non-empty and this allows us to define A(s) = a(s) for s E S. Thus, as a result 
of Proposition 13. ![ we have the following existence result for the modified problem fl3.17p -f l3.2ip . 
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Proposition 3.5. For each ^ G A and each s & S, there is a solution ip = ^/'(l^il, |'C2|,a;3) with 
A = A(|^i|, 1^2!) s) > to the problem with the jump and boundary conditions \3.1^) - ^3.21\) . 

Moreover, ip is smooth when restricted to (—1,0) or (0, 1) with ip{\Ci\, |i^2|,0) 7^ 0. 

Finally, we can use Proposition 13.41 to make a fixed-point argument to find a s G iS satisfying 
s = A(|^i|, 1^2!, s) and to construct solutions to the original problem f l3.2p -( l3^ . 

Proposition 3.6. Let ^ G A, then there exists a unique s G S, so thatX{\^i\, 1^21,-3) = a/— > 
and s = A(|^i|, \^2\,s). 

Proof. We refer to |l8|, Lemma 3.7] for a proof. □ 

Consequently, in view of Propositions 13.51 and 13.61 we conclude the following existence result 
on the problem ([321)-(l33D- 

Theorem 3.1. For each C, E A, there exist ip = ip{\C,i\, |^2|,2;3) and A(|^i|, |^2|) > satisfying 
l\3.^) - (3^^ . Moreover, ip is smooth when restricted to (—1,0) or (0, 1) with ip{\C,i\, |^2|,0) 7^ 0. 

Next, we show some properties of the solutions established in Theorem 13. II in terms of A(^) : = 
A(|^i|, |^2|)- The first property is given in the following proposition which shows that A is a 
bounded, continuous function of ^ G A. To this end, we shall apply the following Ehrling- 
Nirenberg-Gagliardo interpolation inequality, the proof of which can be found in [jj. Chapter 
5]. 

Lemma 3.1. Let Q be a domain in M" satisfying the cone condition. For each eo > there is a 
constant K depending on n, m, p and Eq, such that for 0<e<eo,0<j<m and u G W^'^{Vt), 

V / \D''u{x)\'Pdx <K(e V / |I)"M(x)|Ma; + £-^/(™-^) / 

Proposition 3.7. The function A : A — )■ (0, 00) is continuous and satisfies 

A:=supA(0<7^. (3.31) 



Moreover, 

and if K, > 0, then also 



lim A(n = 0, (3.32) 



lim A(0 = for any e A with \^c\ = (3.33) 
where A denotes the closure of A and \C,\c is defined in Ii3.10\) . 

Proof. The proposition can be shown by a proof procedure similar to that used in 0, Proposi- 
tion 3.9] but with necessary nontrivial modifications in arguments. For the reader's convenience, 
we give the proof in details here. 

(i) We start with proving the boundedness of A. By fl3.24p . we find that 

+ \B?erm" + i^ivi)dx3 - \m9[p\ - «:ien^i(o), 
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which yields 



2/^-iepA(0 £ W^?dx, < \\egWl{Q). (3.34) 



Using the Holder inequality, we can bound 



Jo Jo 



Substitution of the above inequality into f l3.34p gives then 

leP (2/._A(|e|) - Igig]^ l^^pd^s < 0. (3.35) 

Consequently, (I3.35P implies fl3.3ip . since ||'?A^||l2(_i i) > 0. 

(ii) We now turn to the proof of the continuity claim. Since A = a/— a, it suffices to show 
the continuity of a(^) := a(|G|) l^2|)- By virtue of Theorem 13. ![ for each ^ G A there exists a 
function tp^ & A satisfying fl3.2l) - fl3.4p . so that a{^) = E{tp^). Furthermore, ip is smooth when 
restricted to (—1,0) or (0, 1). We have that a{^) < 0, which, when combined with fl3.22p . yields 
the estimate 

I y^/ill^r^? + ^^'rda;3 - l\^\\9[p] - /^I^H^KO) < a(0 < 0. (3.36) 
On the other hand, 

/ \ 1/2 / \ 1/2 A 

mliO) < 2(y^ lerV^ldxa) ( l^^pdxa) < - (3.37) 

because of ips. ^ Thus, plugging f l3.37p into f l3.36p . we get 

I f ^ 5/i||^r^? + rdxs < 2\^\{g[p\ - k\^\')/q-. (3.38) 

Now suppose G A is a sequence so that — )■ ,^ G A. Since A is a open set and ^ ^ (0, 0), 
when n is sufficiently large, there exists a sufficiently small, bounded and open sector domain 
§ C A satisfying the following three conditions: 

• there exists a no > such that ^, ^„ G § for any n > no; 

• (0, 0) ^ S, where § denotes the closure of S; 

• there is a ^ G S such that 

^ < ^, \i\ < lel for any ^ G S- (3.39) 

In order to make use of (I3.38P we have to show that s(^) is bounded uniformly from below 
for n > no. Employing arguments similar to those used in the derivation of fl3.29p . we find that 
there is a G Hq{—1, 1), such that ip ^ and 
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which, together with f l3.39p . imphes 

1 |6^I%|£.|2„7,2 , lJ.n2\A^ 1 



2 



^(1^1'^' + W\')dxs - li9[g] - nl^lWiO) < EoW for any ^ G §. (3.40) 



Let l^lo = inf„>„o |^|, then |^|o > 0. Denoting 



1 



E,{^l^) = -j ^ + 2|el V + 2|^"r)dx3, 

and using (I3.40p . we argue, similarly to f l3.30p . to deduce that 

a{s) < —Ci + SC2 for any ^ G S, 

where the two constants ci = ci(|^|o, Eq^i/j), J(V')) and C2 = C2{Ei{i'), Jdi')) are independent of 
s. Keeping in mind that — = A^(^„) = one gets 

< + C2s(^„) - Ci for any n > hq. 

On the other hand, the fact that ip & A shows that are uniformly bounded in H^{—1, 1) for 
any n > uq. Consequently, by virtue of fl3.39p and (13.180 . the estimate (I3.38P implies the uniform 
boundedness of ip^„ in Hq{—1, 1) for any n > uq. 

Plugging the above boundedness facts on ip^^ into the ODE (13. 5 p in the intervals (—1, 0) and 
(0, 1) respectively, we find that are also uniformly bounded in L^(— 1, 0) and L^(0, 1) for any 
n > riQ. Thus, using Lemma 13.1 [ we infer that ip^^^ are uniformly bounded in if^(— 1,0) and 
if^(0, 1) for any n > tiq. So, up to extraction of a subsequence we conclude that 

strongly in H^oi-^l), ^'(-1,0) and H^{0,1), 

which yields 

a(en)=i?(^cJ^i?(^5)=«(0- (3.41) 
Since (I3.4ip must hold for any such extracted subsequence, one deduces that a(^„) — )■ a(^) for 
the original sequence ^„ as well, and hence a{C,) is continuous. 

Finally, we derive the limits as ^ — )■ when k> 0, and ^ — )■ ^ ^ with = l^lc when n > 0, 
where we restrict the variables ^ in A. By virtue of (I3.24p and (I3.37p . < A^(^) < ^gloW^l/ g_, 
which gives (I3.32p . On the other hand, when k > we may utilize (I3.24p to find that 

which results in (13.331) when k > 0. □ 

Remark 3.2. The stabilizing effect of viscosity and surface tension is evident from the above 
calculations. As shown in [3] , without viscosity or surface tension, there exists a domain D, such 
that A(^) — 7- oo for ^ G © as |^| — ?■ oo. With viscosity but no surface tension, by virtue of 
the definition (13.130 . there are only partial spatial frequencies ^ (including all ^ with = in 
particular) which remain unstable, but the growth of A(^) is bounded. With both viscosity and 
surface tension, only those spatial frequencies ^ belonging to the bounded domain A are unstable, 
and A(^) remains bounded. In addition, in the construction of the normal mode solution to the 
linearized system, the horizontal magnetic field can enhance the maximal growth rate A(^). In 
particular, when = 0, the growth rate reduces to the one for the corresponding equations of 
incompressible viscous fiuids. 
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Remark 3.3. From Proposition 13.71 we immediately infer that 



Ao:=mfA(0>0, 



(3.42) 



where the closed domain D is defined by (13.141) . 



3.3. Construction of a solution to the system li3.2^) - (3^\ ) 

A solution to (I3.5p - (I3.9I) gives rise to a solution of the system (I3.2l) - (l3.4p for the growing 
mode velocity v, as well. 

Theorem 3.2. For each ^ G A, there exists a solution ((^, 6', V^, vr) = ((^(^, X3), 6'(^, X3), '^/'(^jXa), 
7r(^,X3)) with A = A(|^i|, |^2|) > to /i3.^) - l^3^ , and the solution is smooth when restricted to 
(—1,0) or (0,1). Moreover, 



(3.43) 
(3.44) 



Proof. In view of Theorem 13.11 we have a solution (V^, A) := (V^d^il, a^s), A(|^i|, |^2|)) sat- 
isfying f l3.5p -( l3l9|) . Moreover, A > and ip E A is smooth when restricted to (—1,0) or (0, 1). 
Then, multiplying fl3.2P i and fl3.2p 9 by ^1 and ^2 respectively, adding the resulting equations, and 
utilizing (13. 2p ^. we find that vr can be expressed by ip, i.e.. 



vr = vrd^il, 161, X3) = [A/i^'" - (A^^^ + A/i|er + S'eDV'']/(A|er) 
Notice that (I3.2p i can be rewritten as 

^" _ ^x'g + A/x|eP + S'e?)¥'/(A/i) = -^m/^ 
with jump and boundary conditions 

M = 0, MCi^ - = 0, ¥^(-1) = = 0. 

Hence, we can easily construct a unique solution of the form 



ei(cie'^+^« + C2e-»+^^ - Uixs)), on (0, 1), 
ei(c3e"-"^ + C4e-'^-^'3 - /-(xa)), on (-1, 0) 



^ = (6 2:3) = 

to the equation (I3.46P with jump and boundary conditions (13.471) . where 



(3.45) 



(3.46) 



(3.47) 



(3.4J 



a± 



2a±fi± 



TT e 



and 



" Ci " 




1 


1 -1 -1 " 


-1 


C2 










C3 




e"+ 


e-"+ 




C4 







e""- e"- 





)dy. 







/(I) 
/(-I) 
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Similarly to (I33HD, 

^ , f 6(cie'^+^^ + C2e-<'+^3_/^(3,3))^ on (0,1), 

6 := 0[t, X3 = < 3.49 
\ 6(c3e'^-"^+C4e-"-^^-/-(x3)), on(-l,0) 

is a unique solution of f l3.2p 9 with jump and boundary conditions: 

Consequently, {ip,9,ip,Tr) is a solution to the system fl3.2p - fl3.4l) . Now, we define 
((^, e, i), vr) : = (<^(^, X3), 9{^, xs), X3), vf (^, X3)) 

Thus, {ip,9,ip,7!-) is still a solution to the system fl3.2p - fl3.4l) . and moreover, {ip,9,ip,7i) satisfies 

(IS33D. 

Finally, making use of (I3.2p 4 and (13.230 . we conclude that 



1 1 



< J ^(|^r + |^|'+|^|')da;3 



and 

• 1 



-1 



l^'l^dxs < 2/g_ 



The above two inequalities imply (I3.44p immediately. □ 

Remark 3.4. For each 0:3, it is easy to see that the solution (<^(^, ■), 9{C,, ■),ip{C,, ■), 7r(^, ■), A(|,^|)) 
constructed in Theorem 13.21 possesses the following properties: 

(1) A(|,^i|, 1^21)5 ^(^5 ■) and 7r(,^, ■) are even on ^1 or ,^2, when the other variable is fixed; 

(2) (^(^, ■) is odd on ^1, but even on ^2, when the other variable is fixed; 

(3) 9{C,, ■) is even on ^1, but odd on ^2, when the other variable is fixed. 

The next result provides an estimate for the ff'^-norm of the solutions {ip, 9, ip, tt) with ^ 
varying in domain D, which will be useful in the next section when such solutions are integrated 
in a Fourier synthesis. To emphasize the dependence on ^, we write these solutions as (<^(0 = 

Lemma 3.2. Let ^ G D, ^(0 := ^(0; ^iO ■= ^(0; ^iO ■= 7r(0 161) be constructed 

in Theorem \3.^ then for any k > there exit positive constants ak, hk and Ck depending on di, 
d2, Ao, \B\, Q±, fi± andg, so that 



\HH-m + MOWhhoa) + mOWmi-m + 11^(01^^0,1) < a^, (3.50) 

Ih*(-i,o) + llV^(Olk*(o,i) < h, (3.51) 

lk(OllHfc(-i,o) + lk(OlkHo,i) < Cfc, (3.52) 
where di and ^2 are the same constants in 1^3. and Aq is defined by l[3.4^ - Moreover, 

Ml^-i.i) + + = 1- (3.53) 
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Proof. Throughout this proof, we denote by ci, ■ ■ ■ , cs generic positive constants which may 
depend on di, d2, Xq, \B\, g±, fi± and g, but not on |^|. Obviously, fl3.53p follows from fl3.44p 
immediately. 

(i) First we write (13. 5p as 

r'io = [(A'f^ + 2A/iiep + BY,)r{o - \^\\>^'0 + A/xiep + B'ei)m]/M- (3.54) 

If we make use of fl3.3ip . fl3.42p . and the fact |^| < d2, Lemma [3.11 and the Cauchy-Schwarz 
inequality, we see that there exists two constants Ci and C2, such that 

<(c2 + l){s''/^m\\LHi±) + V^\\r'\mLHi,)) forany5G(0,l), 

respectively, where /+ = (0,1) and /„ = (—1,0). Choosing y/e = l/{2(c2 + 1)} in (I3.55P and 
using ( ]3.53p . we arrive at 

iO\\L^ii±) — for some constant C3 > 0, (3.56) 

whence, 

ll^""(0IU^(-i,i)<£4:=2c3. (3.57) 
Thus from (I3.56P and Lemma [3TT] we get 

\\rmLH-i,i)<c,. (3.58) 

and 

||^'"(OIIl^(-i,i)<56. (3.59) 

Summarizing the estimates (13.530 . (I3.44p . and J3.57l) - (l3.59p . we conclude that, for each non- 
negative integer k G [0,4], there exists a constant > depending on d2, \B\, Aq, g±, /i± and g, 
such that 

U^'\mLH-i,i)<h. (3.60) 

Differentiating (13.540 with respect to X3 and using (I3.60p . we find, by induction on k, that 
(Km\ holds for any > 0. This gives (l33T|) . 

(ii) Recalling the expression (13.450 of tt and the fact that |^| > di, we employ (13.420 and 
(I3.3ip to deduce that for any A; > 0, 



which implies (I3.52p . 

(iii) Making use of (K^ . (KW2^ . (ESU) and |^| < ^2, we get 

ll^"(OIU^(-i,i) <^lk(OllL^(-i,i) + (tT^ + ^2 + 11^(0 < cr. (3.61) 

Applying (I3.6ip . (I3.53P and Lemma \3A] we obtain 

WmLH-m + WmiHo,!) <C8. (3.62) 
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Combining fl3.53p with f l3.38p and f l3.62p . we conclude that, for each nonnegative integer k G [0, 2], 
there exists a constant > depending on di, d2, Xq, g±, /i± and g, so that 

y\0\\LH-i,i) < dk. (3.63) 

Thus, by virtue of fl3.46p and induction on k, fl3.63p holds for any A; > 0. Following the same 
procedure as used in estimating ip, we infer that for each k > 0, 

\\0^'\O\\LH-i,i)<dk (3.64) 

for some constant depending on di, d2, Xq, g±, fi± and g. Adding f l3.64p to (13. 63 p . we arrive 
at 

\\^^'\mLH-i,i) + \\0^'\mLH-i,i) < (dk + 4) for any k>0, 
which yields fl3.50p . This completes the proof. □ 

3.4- Fourier synthesis 

In this section we will use the Fourier synthesis to build growing solutions to (ll.20p with 
B = {B,0,0) out of the solutions constructed in the previous section (Theorem 13. 2p for fixed 
spatial frequency C, E A. The constructed solutions will grow in-time in the piecewise Sobolev 
space of order k, &{Q), defined by (12. ip . 

Theorem 3.3. Let B = (5, 0, 0) and f G C^[0, oo) be a real-valued function. For any ^ G M^, 

we define 

wi^, X3) = X3)ei - id{i, 0:3)62 + ^{i, 0:3)63, 

where 

{the solutions 'provided by Theorem \3.2\ 
mt/i A(|ei|,|6l)>o, ^/eero, 
zero, and X{\^,\,\^2\) = Xo, ^/^ ^ ©• 

Denote 



r^{t,x) = — / /(|e|)^(e,X3)e^(l«^l'l«^l)V^'«de, (3.65) 
vit,x) = A /" l6|)/(|el)^(e,X3)e^(l«^l'l«^IV-'«de, (3.66) 



?(^'^) = A / A(|6|,|6l)/(lel)vf(e,X3)6^(l«^l'l«^l)V-'«de (3.67) 

Then, {7],v,q) is a real-valued solution to the linearized problem U.20\) along with jump and 
boundary conditions U.21\) and U.22\) . For every k ^N, we have the estimate 

h(0)||H. + ||t;(0)||H. + ||g(0)||H. <~cJ [ fiim^Y' < 00, (3.68) 

where Ck > is a constant depending on the parameters di, d2, \B\, Xq, q±, h± and g. Moreover, 
for every t > we have r]{t),v{t),q{t) G H^(i7), and 

e*'"h(0)||H. <h(t)||H. <6*^|h(0)||H., (3.69) 

e*^«||t;(0)||H. <||t;(t)||H. <e*^||t;(0)||H., (3.70) 

e*'1|g(0)||H^<||g(t)||H^<6*^||g(0)||H., (3.71) 

where Xq and A are defined in Iji3.42\ ) and Ii3. 31\) . respectively. 
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Proof. For each fixed C. G M^, 



r/(t,x') = /(|el)^(e,a:3)e"(l«^l''«^IV^'^, 

g(t,x) = A(|ei|,|6l)/(lel)vr(e,X3)e^(l«^l'l«^IV-'« 



give a solution to ( ll.20p . Since / G C^[0, oo), by the construction of w, Lemma [3]2] imphes that 

U°°(n) 



sup \\d^w{^, ■)IU°°(n) < oo for all A; G N. 



5Gsupp(/) 

Also, A(^) < A. These bounds show that the Fourier synthesis of the solution given by f l3.65p - 
f l3.67p is also a solution of fll.20p . Because / is real- valued and radial, D is a symmetrical domain 
(see Proposition 13. ip . we can easily verify, recalling Remark 13. 4[ that the Fourier synthesis is 
real-valued. 

The estimate fl3.68p follows from Lemma 13.21 with arbitrary k > and the construction of 
{t], V, q). Finally, we can use 0332]), fl33T]) and ^M^-^Mj to obtain the estimates f l339|) -f l3TT|) . 
□ 



4. Global instability for the linearized problem 

4-1. Uniqueness of the linearized equations 

In this section, we will show the uniqueness of solutions to the linearized problem with lower 
regularity, which will be applied to the proof of Theorem 12.21 in Section 5. For this purpose, we 
need a generalized formula of integrating by parts (or Gauss-Green's formula). Let us first recall 
the boundary trace theorem (see Theorem 5.36 in jll, Chaperter 5]). 

Lemma 4.1. Let U he a domain in M" satisfying the uniform -regularity condition, and 
assume that there exists a simple (m, p)-extension operator E for U. Also assume that mp < n 
and p < q <p* = {n — l)p/ {n — mp). Then, there exists a hounded linear operator 

7^ : W^^P{U) L'^idU), 

such that 

'j^{u) = u on dU 

for allueW'^^''{U)r}C{U). 

The function 7^(m) G L'^{dU) is called the trace of the function of u G iy^'^(f/) on the 
boundary dU. By the Stein extension theorem (see Theorem 5.24 in [l], Chaperter 5]) and the 
definition of the uniform C^-regularity condition (see Definition 4.10 in [l|, Chaperter 5]), it is 
easy to verify that Q, and have different simple (m, p)-extension operators. Keeping these 
facts in mind, we can start to show the following formula of integrating by parts. For convenience 
in the subsequent analysis, we will use the notations 7+(/) := 7^+(/+) and 7-(/) := 7^+(/-)- 

Lemma 4.2. For all v G Ho{^) and w G H-'^(i7), we have 

/ diwvdx = — / wdivdx + / (7+(w) — 7_(w))7+(f )aida; (4.1) 
Jq Jn 

for i = 1,2, 3, where ai = a2 = and ^3 = — 1. 
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Proof. Temporarily suppose v G CliVt), G C^(i7+) and G C^iVt-). By the Gauss-Green 
theorem, we have 

/ diwvdx = — I wdivdx + I — 0)Q:ida;. (4.2) 
JQ Jn JW' 

Using Lemma [4. H one has 

\\{v - -i+{v)){x\ 0)||l2(ir2) <||i; - 7+(^)||L2(af7+) 

= h+{v - t;)||L2(an+) < c\\v - v\\Hi{^a+) 

<c\\'^-v\\Hl{a) 

and 

||(w+ -7+(ii;+))(x',0)||l2(k2) < c\\w+-w+\\H^n+) 
for some constant c > 0. By the Holder inequality, the above two estimates imply that 

-7+(w)7+(t;))(a;',0)||Li(M2) 

< ||(t;(i/;+-7+H))(a;',0)|Ui(M2) + l|(7+H(t)-7+(^)))(2;',0)|Ui(M2) 

< \\v{x\ 0) |U2(M2) II - i+{w)){x\ 0) |U2(K2) (4.3) 
+ ll7+(^)(a;',0)||L2(K2)||(i; -7+(i;))(a;',0)||L2(K2) 

< c^\\v\\Hl(n)\\'^+ - + c^lk+||//i(f7+)||^^ - AhKu)- 
Similarly to (14. 3p . one gets 



(w_i; - 7_(ti;)7+(t;))(a;',0)||Li(M2) 



(4.4) 



In addition, ii Vm ^ v strongly in HliVL), then there exists uiq > such that 

W'^mWuim < WAhK^) + 1 for any m > mo. (4.5) 

Since Co{Q) is dense in Hq{Q) and Co(M'^) dense in H^{Q±), the identity (14. ip follows from 
(I4.2p - (l4.5p and a standard density argument. □ 



Definition 4.1. GivenT > and t/ie initial data (riQ,vo) to the linearized problem U.20\) -in'. 
a triple (ri,v,q) is called a strong solution of M.20\) - 11.22^) . if 

(1) r],ve C0([0,T],L2(fi)), ry(0) = r/o, ^^(0) = vo and 

ess sup (||t;(t)||H3 + + \W)U^ + \Ht)\\min)) < oo. (4.6) 

0<t<T 

(2) The equations 

dtV = V, (4.7) 

gdtV + Vq = fiAv+ J2 BiB„,d^^r], (4.8) 

l<i,m<3 

divi; = (4.9) 
/ioW a. e. m (0, T] x (^] \ {xg = 0}) . 
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(3) For a.e. te (0,T), 

/ (7+(g)/ - ^l+{Vv+ + Vvl) - {^4q)I - fi_{Vv_ + Vv^)))e^ ■ ^pdx' 

= / 9[q]v+3^3 + Y] B^Bi{dir]+ - di7]_) ■ - k V dir]+3di(p3 dx' (4.10) 

-^^^ \ 1=1 r=l J 

holds for any ip G Hq{E?'), where rj^'^ and ^3 are the third component of f]^ and v, 
respectively. 

Remark 4.1. Since v{t) G Hq{Q) fl H^(f2) for each t > 0, we can make use of the embedding 
theorem and (14. 9 p to obtain 

v{t) G c\n), v+{t) G c\n+), v^{t) G c\n^), (4.11) 

v{t) = 0on dn, (4.12) 
V^>v+ = W^'V_ on (4.13) 
divf (t) =0 in for a.e. t > 0. 

Thus, in view of ( I4.13p . we define for the sake of simphcity that 

V^iv := Vx'V+ = Vx'V- on x {0}. 

Moreover, by virtue of Lemma 14. H there is a constant c such that 

||t;(t,x',0)||Hi(M2) < c\\v{t)\\H^n^) for a.e. t > 0. (4.14) 

Remark 4.2. In view the regularity of {ri,v,q) in Definition 14. H we can see that the equahty 
(I4.10p makes sense. 

Remark 4.3. It is easy to verify that any {rj, v, q), which is a solution established in Theorem 
is a strong solution to the linearized system (11. 200 - 01. 220 . 



Theorem 4.1. (Uniqueness) Let B := {Bi, B2, B^) be a constant vector. Assume that {'fi.,v,q) 
and {f],w,q) are two strong solutions to U.2(]\) - fl.2^) . with v{0) = v{0) = vq, 77(0) = f/(0) = rjQ. 
Then {fi,v,Vq) = (r/, Vg). 

Proof. Let {ri,v,q) = {fj — f],v — v,q — q). Recalling Definition 14. 1^ {ri,v,q) is still a strong 
solution to the linearized system (ll.20p -f ri.22p with zero initial data, i.e., r]{0) = and v{0) = 0. 

Let ti G (0,T]. Multiplying (14. 8 P by v, integrating over (0, r) x Q for any r G (0,ti) and 
using (14. Sp . we find that 

/ / gdtv ■ vdxdt + / / dw{ql — fi{Vv + Vv'^)) ■ vdxdt 
Jo Jn Jo Jn , , 

„ r f - - (4-15) 

> / / BiBmOi^T] ■ vdxdt. 
^,^.^nJo Jn 



l<l,m<3 



(1) Firstly, we transform (I4.15P to the form of energy equality. By virtue of the regularity 
dMD, dUD implies that 

dtv G i:2((0,T) X n), 
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which, together with v e L°°{0,T; H\Q)) n C%[0,T], L\Q)), yields 

I I gdtv ■ vdxdt = I- f f)f^(r)da; — ^ f gv'^{0)dx. (4.16) 
Jo Jn ^ Jq 2 Jq 

Thanks to Lemma [4.11 f l4.lip - fl4.12p . and regularity of q, we obtain 

f f div{ql - n{Vv + Vv'^)) -vdxdt 
Jo Jq 

= 11 {l-{q)I - ^l-{Vv^ + Vlr^_) -{^+{q)I - ^l+{Vv+ + Vvl)))e:,■vdx'dt (4.17) 
Jo Jm.^ 

+ / / nVv : {Vv + Vv^)dxdt. 
Jo Jq 

Employing fl4.14p . divf = at the plan {xs = 0} and the arguments similar to those used for 
(14. II) . we can show 



2_] / / nVjdjV±sdx'dt = 0. 
j^iJo Jr'2 



Using the above equality, Lemma 14. H and the regularity of u stated in Remark 14.11 we conclude 
that ^ ^ 

/ / fiVv : Vv^dxdt = 2_, / / P-diVjdjVidxdt = 0. ("4 28) 

Jo Jq i<i,i<3"^o 

In view of fl4.15p - fl4.18p and v{0) = 0, we find the following energy equality 

7: [ gv^{T)dx+ [ [ fiVv.Vvdxdt- f f BiBmdf^r] ■ vdxdt 

2Jn Jo Jq i<i,m<3^o 

- (4.19) 

= / (5'[^']^+3^^3 + B^Bi{diiri+ - dir]_) ■ v - 5i^+35i^3)dx'. 
Jm^ 



Now we continue to transform the above energy equality by replacing 77 with v. Since rj G 
C°([0,T],L2(r])) and r]{0) = 0, the equation (HZD gives 

ri{t,x)= / f (s, a;)ds for any t > 0, (4.20) 
Jo 

which, combined with fl4.13p . yields 

diri{t, x) := diri^{t, x) = diri_{t, x) = / div{s,x)ds, z = 1 or 2. (4-21) 

Jo 

Using fl4.20p . f]4.2ip . Lemma l4T] and the regularity of {ri,v), we deduce that 

/ / '^tmV ■ vdxdt = — / / dmri ■ divdxdt — / {dsr]} ■ vdx'dt 

= - J2 I II d^v{s,x)ds-div{t,x)dxdt (4.22) 
{d^rij ■ vdx'dt, 
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and 



y2 / dir]+sdiV3dx'dt = / / d, 

■^^ Jo JR2 .^^ Jq J^2 Jo 



iVsis, x', 0)dsdiV3{t, x', 0)dx'dt. (4.23) 



Consequently, inserting (I4.2ip - fl4.23p into fl4.19p . we arrive at 

gv'^{T)dx + I I fiVv : Vvdxdt 



Jn 



+ J2 / 

l</,m<3 

9[q\ 




dmv{s, x)ds ■ div(t, x)dxdt 



Jn JO 



Jk^ Jo 



Jr^ Jo 



f3(s, x', 0)dsv3(t, x' , 0)dx'dt 



diVsi^s, x', 0)dsdiV3(t, x', 0)da;'dit 



(4.24) 



(2) Secondly, we further simplify the energy equality to an inequality. With the help of the 
regularity of diVs, the property of absolutely continuous functions and the Fubini theorem, we 
conclude that 



Jq 



9jf3(s, x', 0)ds9jf 3(t, x', 0)dx'dt 



diVs^s, x', 0)dsdiV3(t, x', 0)dMa;' 



R2 ^0 



Jo 
di 



diV3{s,x',0)ds 



dtdx' (4.25) 



diVsi^t, x', 0)dt 



dx' > 0. 



Hence, by we find that 



/ gv'^{T)dx + 2 ^ / / fi\djVi\^dxdt 



JR2 Jo 



Vs{s, x', 0)dsf3(t, x', 0)dx'dt 



(4.26) 



2 ^ BiBm / 




dmv{s, x)ds ■ div(t, x)dx'dt. 



n Jo 



l<l,m<3 

(3) Thirdly, we start to deduce the local- in-time uniqueness from the inequality fl4.26p . Anal- 
ogously to (14.251) . the first integral on the right-hand side of fl4.26p can be bounded as follows. 



'0 Jr^ Jo 

< T 

= 2r 

< T 



IR2 \Jo 



f3(s, x' , 0)vs(t, x' , 0)dx'dsdt 
\v3{t,x',0)\'^dx'dt 

-0 

Usit, x', X3)d3V2,(t, x' , X3)dx3dx'dt 



V3{t,x',0)dt) dx' 



Jm.'^ Ji 



(4.27) 



Jw^ 



rgig] Jo 



\d-Mt)\Mxs + 



rglQ] f 



Jo 



\v3{t)\'^dxs dx'dt 



< 



1 



9[q] 



llVJ^dsV^mhinA 



Jo 
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On the other hand, letting Bq = max{|i?i|, |i?2|, I-B31, 1}, and using the Holder, Minkowski and 
Cauchy-Schwarz inequalities, the second integral on the right-hand side of (14.261) can be estimates 
as follows. 



BiB„ 



l<Z,m<3 




Jn JO 

3 



< 



< 



l<i,m<3 i=l 

«o= E E 

l<i,m<3 i=l 



dmv{s, x)ds ■ div(t, x)dxdt 



dmVi{s, x)ds dx 







n \Jo 

t 



1/2 



\diVi{t, x)\ dx] dt 



1/2 



\Jn 



1/2 

\dmVi{s,x)\^dx ) ds 



1/2 



\div,{t,x)\^dx\ dt (4.28) 



l<i,m<3 i=l 



\Jn 



\dmVi{s,x)\ dx ds 



1/2 



\jn 



1/2 ■ 

\diVi{t,x)\'dx] dt 



Jo, 



\djVi{t, x)\^dxdt. 



Substituting fOTD - fOSD into flCTD . we deduce that 

1<«J<3 ° 

1<M<3 



Now, taking ti = min j^^j, r|, the inequality fl4.29p implies 



\v{t) 



\L^{n) 



< 



vmh(n)dt, TG (0,ti]. 



f^-Q^ Jo 

Moreover, if we apply the Grownwall inequality to f l4.30p . we see that 

ll^Wlli^Cn) = for any r G [0,ti], 
which yields v = 0, i.e., v = v. This, combined with (14.81) and (I4.20p . proves that 

(?7, V, Vg) = {f],v, Vg) for any r G (0, ti] . 



(4.29) 



(4.30) 



(4.31) 



(4) Finally, the local-in-time uniqueness of the solution {fj,v) can be continued onto [0,T]. 
In fact, if ti < T, we can take {'fj,v){ti) = {f],v){ti) as the initial data, and continue the above 
procedure (l)-(3) to obtain the uniqueness of solutions (r/, 5) on (0,^2], where ^2 = 2ti = fil/{2bl) 



or ^2 = T. Hence, after repeating this procedure of extending time interval by finite times, (I4.3ip 
holds for any r G (0, T]. □ 



4-2. Proof of Theorem \2.1\ 
We define 

/3i = c/i + (c?2 - c?i)/3, P2 = di + 2{d2 - di)/3. 
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Fix j>k>0, a>0 and let cj be the constants from Theorem 13.31 For each n E N, let tn 
satisfy 

g2t„Ao ^ ctVc|, (4.32) 

i.e., 

Inc,- ln(Q;n) ^ ^ , / n 
= — ^ + := + Ciln(an), 

Ao 

where Ao is defined by ([332]). Choose /„ G Co"^(M2), such that supp(/n) C B{0, P2)\B{0, Pi), 
where fn is real-valued and radial, and 

fM\n = ^- (4.33) 

Now, we can apply Theorem 13.31 with / = /« to find that [rin(t) , Vn{t) , qn(t)) G solves 
the problem (fT^ - df:^ . It follows thus from (^M) and (lAj) that 022]) holds for all n. 

Recalling supp(/n) C B{(3i,f32) and A(|^i|, |^2|) > Aq, we have, after a straightforward calcu- 
lation and using flX6B|l . f 03|l . flT32|) and (^33]), that 

ll^n(t)|||. > / (l + ler)V*^^l^^l'l^^l)/„2(|^|)||t.(e,a:3)||i.(_i,i)de 



for any t > tn, 



which, together with rjn{t,x) = A(|^i|, \C,2\)vn{t,x), implies (12.31) and (12.41) . This completes the 
proof of Theorem 12.11 

5. Proof of Theorem [m 

In this section we show Theorem 12.21 The main idea of our proof comes from j^, but with 
more complicated computations. We argue by contradiction. Therefore, we suppose that the 
perturbed problem has the global stability of order k for some k > 3. 

Let 6, C2 > and F : [0,5] — ?■ be the constants and function provided by the global 
stability of order k, respectively. Fixing n eN such that n > C2. Applying Theorem 12.11 with 
this n, tn = T/2, k > 3, and a = 2, we find that {fj, v, a) solves (11.201) . satisfying 

||77(0)||h. + ||^(0)||h^. <n-i 

but 



IH3 > 2 for t > T/2. (5.1) 

For £ > we define fj^ = efj^O) and Vq = ^^(O). Then, for e < Sn, ||(^o; ^o)IIh'= < ^- So, 

according to the global stability of order k, there exist 77^, f ^, that solve the perturbed problem 
CT-CT with B = (5,0,0), i.e., 

dti]^ = f^, 

gdtvf + (J,fc - G'^,)dkTi^ = B'dlvt, ^ = 1,2, 3, (5.2) 
divt;^ - tT{GWv^) = 
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with jump conditions across the interface 

Iv^ = 0, mn^j = g[g]4nt + B'{e, + d,^^) ■ n^ld^r^^ + ^^^'^ (5-3) 
and initial data satisfying || (^7o' "^o) llflfe < where 
(GT := (qjLa :=/-(/ + W)"\ 

rf = Ny\N'\ with A^" = (ei + edif) x (ea + ed2f) 

= 63 + €{ei X ^277^ + dif]'' X 62) + e^{dif]'' x (92r7^) (5.4) 
=: 63 + eiV^ 



re 



and 



ei + diT]"^ 


'\'^dlr]^ — 


2(ei + ai77= 


) ■ (e2 + d2if]^)did2r]^ + 62 + <92?7^ 


2^2 


|( 


°i + <9i?7^ 


2 


62 + d2r]'^\ 


p-Kei + ^ir^-) -(62 + 92^)12 



Here we have used the Einstein convention of summing over repeated indices, and to denote 
the both case r^i and 77^ at the interface {X3 = 0} in (15.31) . (15. 4 p and (15. 5p . except for the notation 
\dirjl\ (see Remark [2. ip . Moreover, rj'^.v^.q^ satisfy 

sup II (r/^ ^^ q^) it) ||h3 < F(|| (r/g, v^) Uga) < ^2^11 (r^^ ti^) (0) llga < (5.6) 

0<t<T 

and G C2(f]±) when k > 0. 

Now, defining the rescaled functions fj'^ = rj'^/e, v'^ = v'^/e, (f = q'^/e, and rescahng (15. 6p . one 
gets 

sup \\if,v',q%t)\\^s<l, (5.7) 

0<t<T 

which imphes that 

sup ||(r/^^;^g^)(^)||c;l(n+)+ sup ||(r7^^;^g^)(^)||cl(o_) < ^1, (5.8) 

0<t<T 0<t<T 

for some constant Ki. On the other hand, using the jump conditions (15.30 in term of f^, one has 

sup Wv^mH^.in) < 1- (5.9) 

0<t<T 

Note that by construction, {f]^,v'^){0) = {fi,v){0). Our next goal is to show that the rescaled 
functions converge as e — )■ to the solution {fj, v, a) of the linearized equations (11.200 with 
B = {B, 0, 0) and the corresponding linearized jump conditions. 

5.1. Convergence to the linearized equations 

We may further assume that e is sufficiently small, so that 

sup \\eDf{t)\\Lo.^n) <l/9 and e< 1/(2X2), (5.10) 

0<t<T 

where i^2 > is the best constant in the inequality 

||i^G'||jj2 ^ i^2 II ||h2 II G^||h2 
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for 3 X 3 matrix-valued functions F, G. Then, 

G^={l-iI + eiDff)-')/e 
is well-defined by (IS.lOp and the uniform boundedness in L°°(0, oo; H^(f2)), since 



H2 



n=l 



H2 n=l 



oo oo 



n=l n=l n=l 



(5.11) 



Next, we exploit the boundedness of 77^, v"^, and to control dtrf ^ dtv^ and to give some 
convergence results. The first equation in fl5.2p implies that (9tf/^ = v^, therefore 

sup \\dtfi'{t)U^ = sup ||t;^(t)||H3 < 1. (5.12) 

0<t<T 0<t<T 

By virtue of (15. 7p and (15. lip , the third equation in (15. 2p yields 

lim sup ||divt;"(t)||H2 = lim sup \\eii{G' Dv^){t)\\^2 = {]. (5.13) 

s-^Oo<t<T s-^Oo<t<T 

Expanding the second equation in (15. 2p . one sees that 
gdtvl + d,q' - fiAv^ - B^dl.f, 

= eG),{du{(fh, - Kiljk - eG%)d,vl + (J., - eGl)dkv',)) (5.14) 
- efid,{G'^^djv'^ + G%d,v';) ■= X^ 2 = 1,2, 3, 

whence, by (15.70 and (15. lip . 

lim sup llX^llfii = (5.15) 
e-*-Oo<t<r 

and 

sup ||9(W^||hi < -K^3 for some constant > 0. (5.16) 

0<t<T 

By (15. 7p . (15. 9p . (I5.16P and the sequential weak-*compactness, we see that up to extraction of 
a subsequence (which we still denote using only e), 

{f, v', q') [f], V, q) weakly-* in L°°(0, T; Yi'^in)), (5.17) 

weakly-* in L'^iQ.T- H^{n)) 

and 

dtv' dtv weakly-* in L°°(0, T; Y{\Vt)). (5.18) 
From the lower semi-continuity one gets 

ess sup ||(r7,t;,g)(t)||H3 < 1, ess sup < 1. (5.19) 

0<t<T 0<t<T 

On the other hand, if we use (I5.12p . (I5.16p . and (15. 7p . an abstract version of the Arzela-Ascoli 

theorem (see 1^, Theorem 1.70]), and the Aubin-Lions Theorem (see 12|, Theorem 1.71]), we 
obtain that 

{f]',v') ^ {f],v) strongly in LP (0,T;H2(fi)) and C^{[0,T], L\n)) (5.20) 
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for any 1 < p < oo. This strong convergence, together with the equation dtr]'^ = v^, gives that 

dtV^ dtV strongly in Lp{0, T; lP{Q)). (5.21) 

Rescahng the equations ([El]), utihzing (I513|) -( l515|) . (ISTTj) . (ISTTSj) and dOT]) . the resulting 
equations thus imply that 

dtfj = V, 

gdtv + Va = fxAv + B'^df^f], (5.22) 
divt; = 

holds a.e. in (0,T) x Q. From fl5.20p and the initial conditions {f]'^ ,v'^){0) = {fj,v){0), it follows 
that 

(r],vm = if],vm (5.23) 

as well. 

5.2. Convergence to the linearized jump conditions 

We proceed to show some convergence results for the jump conditions. Multiplying (I5.2p 9 
with ip G (C^((0,T) X rescaling the resulting equation, we find that 

/ / {gdtv' - B^dl.f]') ■ ^dxdt + [ [ {Ijk - eG^^,)vA^dxdt = 0. (5.24) 
Jo Jn Jo Jn 

Integrating by parts and making use of the second jump condition in (15.31) . we conclude that 

/ [ ipidj^dxdt+ [ [ T^d.ipidxdt = - [ [ ipilf^^jdx'dt 
Jo Jn Jo Jn Jo Jk^ 



[ ^mn,jdx'dt+ [ [ ^,{mn,j-m,l)dx'dt 

[ igMvln' + B^ei + edif ) ■ n'ldifj + kH'ti') ■ ^dx'dt 





(5.25) 



+ / / ^m^n,} - {fmx'dt, 

Jo Jr2 

where = H^/e. On the other hand, 

[ [ f^jdjipidxdt = [ [ {fdiYip - n{Vv' + V{v'f) : Vip)dxdt 
Jo Jn Jo Jn 

+ e / ix{G]kdkvl + Glf,dkv])djipidxdt 
Jo Jn 

Putting (I5.24l) -( l5.26l) together, we arrive at 

I I {gdtv' - B^dl^f) ■ ^dxdt + [ [ {fiiVv' + V{v'f) : - fdiv^) dxdt 
Jo Jn Jo Jn 

^,G%dk^dxdt + [ [ f^{G%dkvl + Gl,dkv'^)d,<f4xdt 
Jo Jn 



(5.26) 



T 

Jn 



I I {9[Q]nln' + B'^{ei + edif]') ■ ri'ldif]'} + nH'n') ■ ^dx'dt 

Jo J]R2 

- im)dx'dt. 
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(5.27) 



T 

Jm? 



Next, we deal with the hmit of the equahty f l5.27p as e — )■ 0. Obviously, by f l5.17p . f lS.lSp . the 
bounds (15. 7p and f lS.lip . it is easy to verify that 

lim [ [ {gdtv' - B^df^f) ■ <^dxdt 

+ lim /" [ UiVv' + Viv'f) : Vif - q'diYip) dxdt (5.28) 
Jn 

n{gdtv - B'^dfj^f]) ■ ifdxdt + / (/i(Vt; + Vv'^) : Vv? - qdivip) dxdt, 
A Jo Jn 



and 



lime 

e^O 



f [ ^,G^^,d,^dxdt + r [ ixiq,d,vt + Gld,v^^)d,vAxdt 
Jo Jn Jo Jn 



0. 



(5.29) 



Thus, it remains to analyze the convergence of the two integrals on the right hand of (15.27p . To 
this end, using (15. 4p and (15.50 . we rewrite the first integral as 



/ / igiolVsn' + B\ei + edif) ■ n'ldifj + kH'ti') ■ ^dx'dt 
Jo Jm? 



rT 

lo Jr'^ 




+ 



\e3 + em\ ' {1 + eR')\e3 + em\^ 
^B^F^ + g[QUN^ 



(5.30) 



+ e 



+ 



ifdxdt, 



where 



1=1 



i=l 



- e (9i?72 + 92111 + ^9ir]^ ■ 827]^ 

- 2{dif]2 + + edif]" ■ d2f]^)did2f]l 
+ n' ■N^]{e3 + N^) + A.,,r]IN', 



and 



W = |ei + edifl'd'^f - 2(ei + edif) ■ (es + 682^)8182^ + |e2 + e82r]'\^8lf]' 



Clearly, by the trace theorem, (15.70 and (15.81) in terms of 77^, supo<t<T ||(-R^, -^^)(^)||l°°(m2) 
and suPq<(<^ ||L^()f:)||j^2(]82) are uniformly bounded. Hence, as e — > 0, we have 1 + eR"^ > 0, 
le^ + eN"] > 0, and 



+ 



-> strongly in L°°(0, T; L\W)). (5.31) 



Moreover, by virtue of (I5.20p . (15. 7p and (15. 8p in terms of 77^, using integrating by parts, the trace 
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theorem and dominated convergence theorem, we deduce that 



hm 



{l + eR')\es + em\^ 

2 .T 



dx'dif: 



— > hm 

^ e->0 



JR^ 



[1 + eR')\e3 + em\'^ 



dx'dt 



Ehm 



(Pidi7]ldi 



(5.32) 



dx'dt 



2 .T 

/ / ^i^s^iV^sdx'dt 



=1 







and 



(5.33) 



(5.34) 



In view of (E3ID-(E3SD and fOU]) . we find that 

hm / / (^[^]r/^n" + 5^(61 + edif]') ■ n'ldif]ll + kH'ti') ■ y^dx'dt 

= / (fl'[^']^?3V53 - /«V'5i^?3)5iV53dx'dt. 
FinaUy, using f l5.8p . we easily verify that 

hm f I ^m,n,\ - imWdt = 0. (5.35) 

Consequently, thanks to (15. 28 p . f l5.29p . f l5.34p and fl5.35p . we take to the limit in ( I5.27P to arrive 
at 



T 

Jn 



{gdtv - B^dl^f])ipdxdt + [ f {fi{Vv + Vv'^) : Vv? - gdiv^?) dxdt 

Jo Jq 



T 







(5'[^]^3e3 - K^dir]3diip3)dx'dt. 



(5.36) 



i=l 



On the other hand, similarly to fl4.17p . we multiply f l5.22p 9 by (p and integrate over (0, T) xQ, 
to obtain 



[ [ {Qdtv - B^dl^f]) ■ ifidxdt + [ [ {l^{Vv + Vv'^) : Vif - qdivip) dxdt 
Jo Jn Jo Jq 

= (q+I - fi+{Vv+ + Vvl) - {q^I - fi4Vv- + Vv^))y3-¥>dx'dt, 

Jo J M.^ 



(5.37) 



where we have used the fact that 7±(g) = q± at the interface {xa = 0} since q G H'^(fi) for a.e. 
t G (0,T). Comparing f l5.36p with f l5.37p . we conclude that 

/ (q+I - i2+{\/v+ + \/vl) - {q-I - /i_(Vt;_ + Vt;^)))e3 ■ 0dx' 

= / (^[^fe03-«;$^5ir/35*03)dx' for a.e. t e (0,T) (5.38) 
holds for any G H^{M.'^). 
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5.3. Contradiction argument 

Notice that f}^ in fl5.38p includes both cases of 17+3 and f]_^. In view of Definition I4.H we 
find that {f],v,q) is just a strong solution of the linearized problem f ll.20p -f ll.22p . By Remark 
14.31 {fi,v,q) constructed in Theorem 12.11 is also a strong solution of fll.20p - p.22p . Moreover, 
77(0) = 77(0) and v{0) = v{0) (see f l5.23p ). Then, according to Theorem 14.11 

v = v on [0,T) X n. 

Hence, we may chain together the inequalities fl5.19p and fl5.ip to get 

2 < sup ||^'(i)||H3 ^ sup \\v\\ff^3 < 1, 

T/2<t<T 0<t<T 

which is a contraction. Therefore, the perturbed problem does not have the global stability of 
order k for any A; > 3. This completes the proof of Theorem 12. 2[ 
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